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Abstract 



We derive and evaluate expressions for the dc tunneling conductance be- 
tween interacting two-dimensional electron systems at non-zero temperature. 
The possibility of using the dependence of the tunneling conductance on 
voltage and temperature to determine the temperature-dependent electron- 
electron scattering rate at the Fermi energy is discussed. The finite electronic 
lifetime produced by electron-electron interactions is calculated as a func- 
tion of temperature for quasiparticles near the Fermi circle. Vertex correc- 
tions to the random phase approximation substantially increase the electronic 
scattering rate. Our results are in an excellent quantitative agreement with 
experiment. 
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I. INTRODUCTION 



The development of high mobility double-quantum-well structures and techniques to 
independently contact the two wells even when separated by only tens of nanometers, have 
together opened a new fruitful area of research in low-dimensional physics. For example, 
inter-layer electron-electron interactions have been investigated through the frictional drag 
voltage occurring when charge in one layer is moved relative to charge in the nearby layerBi. 
Separately contacted double-well structures have turned out to be extremely useful not only 
in studies of properties unique to bilayer systems but also for examining some properties of 
a single-layer 2D electron gas. In the case of a gated heterostructure with remotely spaced 
quantum wells, where the inter-layer interactions could be safely neglected, Eisenstein et 
a/.i were able to relate the electric field leaking between layers in response to the voltage 
change on a remote gate to the compressibility of the electron layer closer to the gate. 
2D-2D tunneling!'! provides another example where double-well systems can be used to 
probe the electronic properties of individual electron layers in new ways. An important 
feature of ideal 2D-2D tunneling is the conservation of electron momentum; for 3D-3D 
tunneling the component of momentum in the direction of the tunneling barrier is not 
conserved. Momentum conservation is not perfect because of disorder in the tunneling 
barrier and in the electron layers, and also because of inelastic electron-phonon and electron- 
electron scattering. An electron conserving its energy and momentum can tunnel only when 
the subband edges of the two layers are precisely aligned, resulting in large peak-to valley 
ratios in the observed resonant tunneling peaks. Experimental studies in which the Fermi 
surfaces of the 2D electron systems were mapped by measuring the tunneling conductance 
in magnetic fields applied parallel to the 2D planesll have provided, arguably, the most 
striking demonstration that electron momentum is conserved to a remarkable degree in 
GaAs/AlGaAs double- well structures. The heights of resonant peaks in 2D-2D tunneling 
conductances are limited by scattering processes which do not conserve the momenta of 
individual electrons. At zero temperature, peak heights and widths may be used to measure 
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the elastic scattering rate due to disorder.! It has been suggested!!! that the temperature 
dependence of 2D-2D tunneling conductance peak heights and widths can be used to probe 
inelastic electron-electron or electron-phonon scattering processes in the individual 2D layers. 
By comparing with the measured mobility of their system, Murphy et a/.! concluded that 
electron-phonon scattering could not account for their observations and suggested that the 
experiment provides information predominantly about electron-electron scattering rates. 

In this article we discuss the connection between electron-electron scattering rates in 
individual layers and 2D-2D tunneling.! We find in Section II that the energy dependence of 
the electron-electron scattering rate makes the relationship to tunneling conductance more 
complicated than in the case of elastic disorder scattering. The formal expressions derived 
in Section II nevertheless allow the tunneling conductance to be calculated from the energy- 
dependent quasiparticle lifetime. The problem of calculating the quasiparticle lifetime due 
to Coulomb interactions is a standard one in the many-body theory of the electron gas; for 
most purposes the random-phase-approximation (RPA) is reasonably accurate. For the two- 
dimensional electron gas, however, confusing disagreements exist among various analytic 
evaluations of the RPA expressions.0~0 To clarify the situation, we present in Section 
III a detailed derivation of the approximate analytic formulas for the temperature and 
energy dependence of the electron-electron scattering rate. We follow a line similar to that 
described, in detail, in Ref. |ll| and emphasize the points where present work departs from 
previous studies. In Section III, we also discuss calculations of the quasiparticle lifetime 
which go beyond the RPA by including local-field corrections to the effective electron-electron 
interaction. The corrections approximately account for density and spin-density correlations 
present in the ground state of the interacting electron system. When they are included in 
lifetime calculations along with the energy-dependence of the electron-electron scattering 
rate , excellent agreement is obtained, as we discuss in Section IV, with the 2D-2D tunneling 
experiments of Murphy et al3 In Section V, we briefly summarize our results. An account 
of a preliminary version of this work has been presented previously.0 
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II. 2D-2D TUNNELING CONDUCTANCE FOR T / 



We consider a GaAs/AlGaAs heterostructure with two identical quantum wells and equal 
layer densities. For typical separations between the 2D-layers (~ 300A) inter-layer interac- 
tion effects, including the screening of the intra-layer Coulomb potential by electrons from 
the opposite layer, are weak. All such effects will be ignored in following calculations. The 
Hamiltonian can be written as the sum of three terms: 

H = Hr + Hl + Ht 

Ht = - E + h - c -) • (1) 

Hr and Hl are the Hamiltonian for isolated electrons in right and left wells including, in 
general, contributions from intra-layer interactions and from disorder in each layer. The 
tunneling Hamiltonian, Hj<, couples the two systems; we assume that this term can be 
treated at leading order in perturbation theory. For tunneling barriers which are invariant 
under translations perpendicular to the barrier tt p is zero for k ^ k' and is independent of 
k; t^p = tS^p. For noninteracting electrons t determines the difference in energy between 
symmetric and antisymmetric combinations of subband states for the two-layers. 

The Kubo formula, which we will use to calculate the tunneling current /, treats the 
Ht to leading order in perturbation theorylll For 2D-2D tunneling our theory is valid to 
leading order in 2tr/nfi where r is the lifetime of an electron in the individual layers. For 
noninteracting 2D-electrons in disorder free double-quantum-well systems it is never valid 
to treat Ht as a perturbation. The condition for the validity of the weak perturbation 
assumption is that the mean time for an electron to move between layers (r t = irh/2t), is 
much longer than the lifetime of electrons due to scattering within a layer, r. An electron 
hopping from one well to the other will then scatter, i.e. change its in-plane momentum, 
many times before it jumps back to the first well. 

Following a familiar 

line§@ 

of reasoning we obtain 
I(V) = -^t 2 S J ™J_ B ^A(E,k)A(E + eV, k) [n F (E) - n F (E + eV)} , (2) 
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where S is the area of the two-dimensional electron systems, eV is the difference between 
chemical potentials in the right and left quantum wells, rip(E) is the Fermi distribution 
function and A(E, k) is the spectral function related to the retarded Green's function by 

A(E, k) = -2lmG ret (E, k) = -2ImS^(E, fc) ^ 

E - £ fc - ReE ret (E, k) + lm£ ret (E, k) 

We choose to measure energies from the Fermi energy so that = h 2 k 2 /2m — Ep. Near the 
quasiparticle peak the spectral function can be approximated by a Lorentzian: 

Here we have neglected the real part of the self-energy which leads to a physically unimpor- 
tant rigid shift in the quasiparticle energies, causes the quasiparticle effective mass to differ 
slightly from its free electron value, and slightly reduces the weight of the quasiparticle peak 
in the spectral function. These effects play a minor role in 2D-2D tunneling experiments^ 
and we neglect them here in favor of the main effect which comes from the broadening of 
the quasiparticle pole. The width of the Lorentzian peak is related to the self-energy by 
r(6fc, k) = -2Im£ ret (£ fc , k). 

In the limit of r(^, k) — ► 0, corresponding to the noninteracting disorder- free 2D-electron 
gas, Eq.(^) can be written as 

Op f)„ POO rlrp rCO 

I(V) = -t 2 S^ / — 5(x)8(x + eV) / d£ h [n F (x + £ fc ) - n F (x + £ k + eV)] , (5) 

where go — rn/fth 2 is the free-particle density of states for a 2D electron system and x = 
E — Because of the 5-functions in Eq.(|^) the integral over ^ gives eV and we obtain for 
\eV\ -C Ep that the tunneling conductance G(V) = I(V)/V is proportional to S(V). This 
sharp voltage-dependence of the tunneling conductance is a direct consequence of electron 
energy and momentum conservation during the tunneling. An electron with kinetic energy 
£/c can tunnel only when its potential energy is conserved, i.e., when the energy levels in the 
quantum wells are aligned. For identical wells this condition is satisfied only at zero voltage. 
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Scattering processes lead to an uncertainty in the energy of an electron with a given 
momentum k resulting in a finite width of the spectral function and broadened peaks in 
the G(V) curve. At zero temperature the broadening is dominated by elastic scattering. 
Replacing the spectral width r(^, k) in Eq.(^) by h/r e i, where r e i is the constant elastic 
scattering lifetime, Eqs. (0) and (|]) give for h/r el <C E F and \eV\ <C E F : 

° {V) ~ h l S Y(eVy + (h/r el r ■ (6) 
We define Tq to be the half width of the G(V) curve. The dependence of the tunneling pa- 
rameter t on the bias voltage can, and will, be neglected. To see this note that t ~ exp(— kcLb) 
where ds is the width of the barrier between the quantum wells, n = (2mVB/fr 2 ) 1//2 gives 
the decay rate of the wavefunction in the barrier, and Vb is the barrier height. When a bias 
potential is applied the average value of Vb is changed by ~ eV. Using dn/dVB ~ k/Vb the 
t should change by a factor of ~ exp(— ndB{eV) /Vb)- In the experiments to which we refer 
Vb ~ 300 meV and the maximum value of \eV\ is ~ O.lEp ~ 0.5meV. We see that the mag- 
nitude of the argument of the exponential above is much smaller than one. It follows from 
Eq. (Q) that Tq = h/r e i for a system with only elastic scattering. It has been established 
experiment allyili that Tg is temperature dependent, indicating that some inelastic scatter- 
ing process is contributing to the quasiparticle scattering rate. Samples with different levels 
of disorder have G(V) peaks whose half-widths appear to differ by temperature-independent 
constants. This property enables elastic and inelastic contributions to the scattering rate 
to be separated experimentally; all the calculations in this article are for a disorder-free 
system. Our calculations will help confirm the experimental analysis of Murphy et al. who 
attributed the temperature dependent broadening of the peak in the G — V characteristic to 
electron-electron interactions. If we approximate the spectral width due to electron-electron 
interactions, r e _ e (£fc, T), by r e _ e (0,T) we arrive again at Eq.(||) with h/r e i replaced by 
r e _ e (0,T). We show later that this neglect of the energy dependence of r e _ e (£ fc ,T) under- 
estimates T G . 
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III. ELECTRON-ELECTRON SCATTERING RATE 



As explained above, both the temperature dependence and the energy dependence of 
r e _ e need to be calculated in order to compare theory to measured G — V characteristics. 
In this section we present calculations for a pure 2D-electron gas in the random phase 
approximation and in the local-field-corrected RPA. The contribution of electron-electron 
scattering to the spectral width of the one-particle Greens function can be expressed!^]'! in 
terms of the scattering rates of electrons (h/r e ) and holes (h/r h ): 

r e _ e (£ fc ,r) = — — - + 



T e (£ k ,T) T h (£ k} T) 

h _ r d 2 k' r d 2 P 



T h (£ k ,T) [l-n F (^)]r e (^,T) " U 
In Eq.([7p, k,p are the initial electron momenta and k',p' are the final electron momenta. 
Because of the conservation of the total momentum in electron-electron scattering processes, 
p' = k + p — k' . W aa ' is the scattering function which we now discuss. 



A. Scattering function in the random phase approximation 



In the RPA, the scattering function W aa> is spin-independent (W PPA = W PPA = Wrpa) 
and depends only on the momentum transfer q = \k' — k\ = \p' — p\ and energy transfer 

hu: = £ k ' — 6c = £p — Cp' : 

v(q) 



W. 



RPA 



2tt 



(8) 



£rpa(<?, M 

In Eq.(||) v(q) = g^qrFll is the unscreened Coulomb interaction and £rpa(<2S tw) = 
1 — v(q)xo(q,fi<u) is the RPA dielectric functionS0'0 for the 2D-electron gas. (Here 
Qtf = 5'oe 2 /2e is Thomas-Fermi screening wavevector and xoilj^) is the susceptibility 
of a noninteracting electron gas@l^Hlil.) 



At low temperatures the energy transferred in electron-electron scattering (~ kpT) 
is small compared to the Fermi energy and the magnitude of the momentum transfer 
is restricted to the interval (0,2kp). In this limit the real part of the susceptibility 
Rexo(q, hoj) ~ —go, the imaginary part lmxo(q, hu) « 0, and we can write 

W RPA « — tt 1 - — • 9 

Previous analytic evaluations of the RPA quasiparticle lifetimes have employed the approx- 
imation Wrpa ~ 2-n j g 2 , which would be reasonable if qpF 2kp or if the slowly varying 
function Wrpa{q) is multiplied, in the integrand in Eq.(^), by a function sharply peaked near 
q = 0. Neither of these assumptions is valid however since: (i) The condition qpF 3> 2kp 
corresponds to the dimensionless parameter r s = (qTF /2kp) / \/2, conventionally used to 
render the density of an electron gas, being much larger than 1. At such low electron den- 
sities it is known that the RPA fails. In fact, in the next section it is shown that even for 
r s ~ 1 vertex corrections to the RPA substantially increase the electron electron scattering 
rate. Moreover, in GaAs, for electron concentrations typical of tunneling experiments the 
Thomas-Fermi screening wavevector is comparable to the Fermi wavevector; (ii) As we now 
discuss in detail the integrand in Eq.(|7|) has sharp peaks near both forward (q = 0) and 
backward (q = 2kp) scattering momentum transfers. 

In polar coordinates d 2 k' = k' dk' d9y and d 2 p = pdpd8 p , where the angles Qy and 9 p are 
measured with respect to momentum k. The integral over 9 P may be performed using the 
energy-conservation 5-function. Since £ p / depends on 9 P we may write 



09 p 



S{9 P - 9 P ,) , (10) 



where 9 P ^, is an angle at which both the energy conservation 

p' 2 = e + p 2 - k' 2 (ii) 

and momentum conservation 

p' 2 = (k+p — k') 2 = k 2 +p 2 + k' 2 + 2kp cos 9 P - 2kk' cos 9 k > - 2pk' cos(9 k > - 9 P ) (12) 



conditions are satisfied. From Eq.(^) it follows that 
dp' 2 



de r , 



—2kp sin 9 p + 2pk' sin(Q p — 6k') = —2pk' sin 9i 



cos (6>p + z) 

r 

COS z 



(13) 



where 



tanz 



k — k' COS Qui 



k' sin Qfri 

The combination of Eqs. ([TT|) , (|l^) and QUI) gives 



(14) 



2 pk'smQ k> . , 
/e = sml^p + zj + /e/e cosQk< 



cosz 

and together with Eq . flT3|) we finally obtain 

-l 



(15) 



dQ r , 



A+{E F + £ k ){E F + £ k ,)sm 2 Q k 



-1/2 



A = (f P - &<) [6 + 6' + 2£f - (6 + ^f) 1/2 (^' + Ep) 1 ' 2 cos fc , 



(16) 



In the limit of small T and the Fermi functions restrict energies (measured from the 
Fermi energy) of particles involved in the scattering process to a small region near zero 
energy. As seen from Eq. ([l6|) , A is then small and the integrand in (0) has equivalent sharp 
peaks in the available phase space for scattering near Qy = and Q k > — tt. The main contri- 
butions to the electron-electron scattering rate come from processes with small wavevector 
transfer (forward scattering) and wavevector transfer ~ 2k p (backward scattering). This 
suggests that Wrpa^) m &y be approximated by the average of its value at q = and its 
value at q — 2kp\ 

2tt 
7o 



W RPA « —w f p b PA , 



(17) 



where 



w 



f,b 
RPA 



l+fl + ^ 



v/2 



(18) 
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B. Scattering function in the local-field-corrected RPA 



In estimating quasiparticle scattering rates the RPA does not account for electronic 
correlations in the interacting electron gas. Technically, the RPA for the electronic self- 
energy neglects vertex corrections to the dynamically screened exchange energy. On physical 
grounds correlations are expected to suppress scattering between like-spin electrons since 
these electrons are required to avoid each other by the Pauli exclusion principle and to 
enhance scattering between opposite-spin electrons. In general the four-point scattering 
amplitude for electrons depends on the energies and momenta of all electrons and not just 
on the momentum (q) and energy (u) transferred in the scattering event. Nevertheless, 
a number of workers^!! have suggested similar approximations in which the four-point 
scattering wavefunction is replaced by an effective electron-electron interaction dependent 
only on q and uj. In these approximations correlations are accounted for by modifying the 
bare electron-electron interaction to take account of the correlation clouds carried around 
by each quasiparticle. 

Here we estimate corrections to the RPA by adopting the effective electron-electron 
interaction suggested by MacDonald and Geldart.0 Their effective interaction has density- 
density {t nn ) and spin-spin (t mm ) components: 

^LFRPA = ^{tnn + tmm) 

^LFRPA = 27r(t nn — t mm ) . (19) 

Here the density-density interaction, 

, = v{g)+F nn (q) 

nn l- X o(q,huj)(v(q) + F nn (q)) ' 1 J 

is a screened Coulomb interaction and the spin-spin interaction, 

j. F mm (q) 

can be thought of as exchange interaction which favors parallel spin alignment. The local 
field factors F nn (q) and F mm (q) which appear in Eq.flSOp and Eq.(pT|) are related to the static 
density (Xnn(q)) and spin (x mm (?)) response functions of the electron gas: 
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Fnn(q) = Xo(q) 1 - Xnn(q) 1 
Fmm{q) = Xoiq)' 1 - Xmm{q)~ l - (22) 

From Eq.(|2*2]) we see explicitly that the local fields vanish in the RPA. Quantum Monte- 
Carlo calculations of the response functions indicate that the wavevector dependence of the 
local fields is weakil for q < 2k p. Here we approximate the local fields by their q — > limits 
which are@ accurately known. Using the same arguments as in the previous section, the 
four point scattering function W[p RPA can then be approximated by the average between 
its forward and backward scattering limits, i.e, 



Wl FRPA + W LFRPA _ 2 2 2tt Lb 

— Z7V \ l nn + l mm) ~ o LFRPA > \ ZO ) 

9o 



where 

n J,b _ 1 + ( 1 + r 3 V2+g F nn ) ( ggjkm \ (r)A s 

WLFRPA ~ 2 + {l + 9 oF mm ) ■ (24) 

In Figure |I| we compare scattering functions in the random phase approximation and 
in the local-field-corrected RPA as a function of density. The RPA result is recovered in 
the high-density limit where the local fields go to zero. For r s = 1, corresponding to the 
electron density of the sample in the experiment^ of Murphy et al, the local- field-corrections 
increase the scattering rate by approximately 30%. As also shown in Figure [I], there is a 
~ 50% difference at r s = 1 between the forward scattering approximation for the scattering 
function and the approximation obtained by averaging forward and backward limits. (In 
the forward scattering limit w f RPA = 1 and w{ FRPA = 1 + [g F mm /(l + g,F mm )}\) 

C. Approximate analytic results for low temperatures and energies 

Analytic expressions for the RPA electron lifetime h/r e at low energies and/or temper- 
atures have been derived previously in four independent studies^FEl with slightly different 
results obtained in each case. None of these expressions agree with the analytic results we 
present below, whose accuracy has been verified by comparing with independent numerical 
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calculations. While some of the discrepancies appear to be due to inadvertent algebraic 
errors, some are associated with significant aspects of the physics of electron-electron scat- 
tering in two-dimensional electron systems which we draw attention to below. Because of 
the existing confusion we give a detailed description of our calculation. We have found the 
detailed analysis presented in Ref. [TT| to be very helpful and have followed this calculation 
closely. 

We have shown above that the phase space for electron-electron scattering at low- 
temperatures is dominated by contributions with equal weight near the forward and back- 
ward scattering limits. This property invalidates the approximation, made in previous an- 
alytic evaluations, in which the RPA screened interaction is approximated by its forward 
scattering limit. For our analytic calculation we have adopted an approximation which is in 
the same spirit by replacing the scattering function by Eqs.([T7|),([T8|) or Eqs.(p^),(p4[). Mak- 
ing this replacement we can take the interaction strength outside all integrals. We integrate 
first over the angle Qy between the incoming and scattered momenta in Eq. (|7|). Using 
equations similar to (|T^-(T^) we find that 



- Qitk - 60 + (E F + £k){E F + Q sin 2 9 k ) 1/2 . (25) 



de v 

As illustrated in Figure |2| there are two different angles, 9 k ',i and dk',2, with identical inte- 
gration weights ( PS| ) which satisfy the energy and momentum conservation conditions. The 
corresponding momenta are related by a mirror symmetry with respect to the vector k + p. 
(The associated factor of two in the scattering rate appears to have been missed in some 
previous work.) At low T we can take 

(E F + £ k )(E F + QnE F . (26) 

Defining «f fc = i k jk B T, & = £ k >/k B T, and £~ = i v jk B T leads to 
h m 2 (k B T) 2 2vr -- f°° 



/OO _ /"CO _ _ 

-oo J —CO 



r e (a,T) (27r/^ E F g 2 

1 - n F (( v )] [l - n F (( k + i P - eXOl 1(B) , (27) 
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where w^ h stands for either w^ PA or w^ b FRPA and 



B 



E P 



(28) 



The angular integral 1(B) reads 



1(B) = 2 



7T-0 O 



d6 r , 



00 2{B + sin 2 0, 
du 



1/2 



«o (1-M 2 ) 1/2 (5 + M 2 ) 1/2 



(29) 



Here uq = for i? > and wo = y for B < 0. To find the leading low-temperature 
behavior of the electron-electron scattering rate we replace the expression (B + u 2 ) 1 ' 2 by 
first two terms of its Taylor expansion and approximate the elliptic integral Eq.(^9|) by 



1(B) » ln(8) - In IS) 



(30) 



Inserting Eq.(|3"0D into Eq. (B7[) and using the identity 



d£u I di P n F (Q \l - n F (( k ,)] \l - n F (i k + £ p - („)] = \(-k 2 + £ fe 2 ) [l - /(&)! (31) 



we arrive at 



where 



h/r e (£ k ,T) wf' b (k B T 



E f 



Ei 



r 2, 



1 - ni?(Cfc 
In 8 , ( E F 



- ^(6 



(32) 



^(4 



1 r 



1 - n F (i k ) 



— 1 Z" 00 



k' I ""Sp - 
oo J — oo 



x n F (Q l-n F (Ck') l-n F (i k + L-i 



(33) 



In the limit of — > 0, i.e., for an electron on the Fermi surface, F(0) = 0.41388 and we 
obtain the leading contribution to the electron-electron scattering rate at low temperatures 



h/r e (0,T) 
E F 



fb n ( k B T" 

« w 1 ' — 

4 I Ep . 



E F \ In 8 

In — - + .083 

k B T 2 



(34) 
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As follows from Eq.flT]) the scattering rates of electrons and holes are equivalent at = 0, 
i.e., 

n 



r e _ e (0,T) = 2- 



(35) 



Te(0,T) ' 

At zero temperature the calculation is simplified by the fact that the Fermi distribution 
functions are reduced to step functions. We find that for > 



Ep Ep 
and for < 

r e _ e fe,o) tyr fc (&,0) 



2vr V £ F 



In 



w 



2tt 1 E F 



In 



.6 



Ep 



In 8 

+ + -5 



In 8 „ 
+ - + .5 



(36) 



(37) 



IV. NUMERICAL RESULTS 

All the calculations we discuss in this section have been performed for r s — 1 . In 
Figure |3] the low-temperature and low-energy approximate analytic results (Eqs.(|3"4"D,(55) 



and (P^) , (P?D ) are compared with numerical results obtained within the constant interaction 
approximation (|2"3"|) , (|2~4]) . For temperatures or excitation energies up to 10% of the Fermi 
temperature or Fermi energy, resp., the terms with higher powers of T or ^ can be safely 
neglected. The discrepancies here are due only to the low-temperature and low-energy 
approximations. We see that including the T 2 and £| corrections to the leading — T 2 ln(T) 
and — £jiln|£fc| terms greatly extends the range of validity of these expressions. Note that 
the spectral width due to electron-electron interactions is not precisely an even function 
of energy. In the constant interaction approximation at zero temperature the difference 
between scattering rates of an electron with energy > and of a hole with energy — is 
due to the term given by Eq. (|25|) . Since (Ep + |^|) > (Ep — |^|) the integrand in Eq.(^) 
is smaller for electrons than for holes with the same absolute value of energy. Making the 
approximation of Eq.QSBj) which only affects terms of higher order in than £| In |^| and 
£f gives the leading low-energy behavior which is even in 
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The validity of the constant scattering amplitude approximation (Eq.(|23"D and Eq.(|24l)) 
was examined by performing three different numerical calculations of r e _ e (0,T) and 
r e _ e (£fc, 0). In Figure £| we show results obtained using (i) the constant scattering function 
W[ FRPA = 27r/gQW^p RPA (These results are identical to the curves labeled "Numerical" in 
Figure [|), (ii) the g-dependent approximation W a (q) to the scattering function which results 
when the wavevector and frequency- dependence of the susceptibility Xo(<L is neglected 
(Xo(l,u) ~ — go), and (hi) the full q and uj dependent scattering function (WiFRPAiq, ^)) in 
the local-field-corrected RPA where both the real and imaginary parts of xoiq,^) are taken 
into account. We see that averaging between forward and backward scattering limits pro- 
vides an excellent approximation for the (/-dependence of W a (q). The discrepancy between 
r e _ e calculated with the constant iW[ FRPA ) and full scattering function (W[,FRPA{q, w)) is 
primarily due to neglecting the imaginary part of xo which is non-zero at finite frequencies. 

At T ^ quasielectrons and quasiholes are both present both above and below the Fermi 
energy. The energy- dependent spectral width r e _ e (^) is then a sum of electron and hole 
scattering rates. Numerical results showing electron and hole scattering rate contributions 
are shown for T = 0.3 Tp, in Figure Spectral widths calculated at several different 
temperatures are plotted in Figure || Dotted curves show the results of calculations where 
the temperature dependence of the susceptibility xo is neglected. This approximation greatly 
simplifies non-zero temperature numerical calculations, however, we see here that it can be 
safely used only for temperatures T < 0.1 T F . 

Inserting r e „ e (^.) into Eq.(|4]) we can evaluate Eq.(|^) for the tunneling current and cal- 
culate the tunneling conductance as a function of applied voltage at a fixed temperature. 
At low temperatures and energies the spectral width T e _ e (^k,T) is, roughly, the sum of 
r e _ e (0,T) and r e _ e (£ fc ,0) (see Figure ^|). Looking at the approximate analytical formulas 
(@)-(0) we see that if T/T F < 1 and \£ h \/E F < 1 then also T e _ e (£ k ,T)/ E F < 1. For 
\eV\/E F < 1 we can therefore replace the difference between Fermi functions in Eq.(fJ) by 
dn F {E) / BE x eV resulting in the following formula for the tunneling conductance 
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ay) = 2 ^s 9 4 r ™ r ,i e -f k \ w , a 



ft 2 J-oo 2tt (£ - £ k y + (r e _ e (£ fe )/2)< 

T e ^ k -eV) dn F (E) 



(38) 



(E + eV — a) 2 + (r e - e (6 - eU)/2) 2 <9£ 
Using Eq.(p8|) we determined numerically the effect of the energy dependence of the 
spectral width due to electron-electron interactions on Ta(T). As shown in Figure |7|, ap- 
proximating r e _ e (£fc,T) by r e _ e (0,T) underestimates Tq(T) by a factor ~ 1.1-1.16 in the 
range of temperatures used in the tunneling experiment!. In Figure || the calculated Tg{T) 
is compared to the measured half-width of the tunneling conductance peak, showing ex- 
cellent agreement between the theoretical result obtained in the local-field-corrected RPA 
and experimenal data. Note that the agreement between theory and experiment has been 
achieved without introducing any adjustable parameters. As for the scattering rate calcula- 
tion, local-field-corrections increase the RPA Tg(T). At T jTp = 0.15 the correction factor 
is ~ 1.3 and improves agreement between theory and experiment. While it is certainly 
possible that the excellent agreement with experiment is partly fortuitous, our results leave 
little doubt that the temperature dependence observed in the 2D-2D tunneling experiments 
is due to the finite quasiparticle lifetimes which result from electron-electron scattering. 



V. SUMMARY 

In 2D-2D tunneling experiments, resonant tunneling features have no width in the ab- 
sence of quasiparticle scattering processes. This property of 2D phase space allows 2D-2D 
tunneling experiments to measure quasiparticle scattering rates. In this article we have 
examined the role of electron-electron scattering in recent 2D-2D tunneling experiments. 
For elastic scattering the width of the resonant tunneling peak is proportional to the scat- 
tering rate for electrons at the Fermi energy. For electron-electron scattering we find that 
the relationship is complicated by the non-negligible dependence of the scattering rate on 
quasiparticle energy relative to the Fermi energy. We have shown that, although the energy 
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dependence of the scattering rate results in a non-Lorentzian lineshape for the resonance, it 
has only a modest influence on the temperature dependence of the half-width of resonant 
tunneling peak. Analytic expressions have been derived which approximate RPA and local- 
field-corrected RPA results for the scattering rates at low temperatures and energies. The 
accuracy of these analytic expressions has been confirmed by comparison with independent 
exact numerical evaluations of the scattering rates in these approximations. For the electron 
densities of existing experiments local-field-corrections increase the RPA scattering rate by 
a factor ~ 1.3 and are important in improving agreement between theory and experiment. 
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FIGURES 

FIG. 1. Forward (/) and the average between forward and backward (/, b) limits of the scatter- 
ing function calculated in the random phase approximation (RPA) and in the local-field-corrected 
RPA (LFRPA). 

FIG. 2. Schematic illustration of angles between the incoming and scattered momenta for 
two-dimensional electron-electron scattering. 

FIG. 3. Spectral width as a function of temperature at zero energy (upper panel) and as 
a function of energy at zero temperature (lower panel) calculated using the constant scattering 
amplitude approximations discussed in the text. These results are for the local- field-corrected 
RPA. Numerical calculations (full lines) are compared to these approximate analytic results. 

FIG. 4. Spectral width as a function of temperature at zero energy (upper panel) and as a func- 
tion of energy at zero temperature (lower panel) calculated numerically in the local-field-corrected 
RPA. Exact LFRPA calculations (full lines) are compared to results obtained with (/-dependent 
approximate scattering function (dashed lines) and with the constant scattering amplitude (dotted 
lines) . 

FIG. 5. Spectral width (full line) as a sum of electron and hole scattering rates (dotted lines) 
at T = 0.3 Tp. 

FIG. 6. Spectral widths (full lines) calculated at different temperatures. Dotted lines show 
results obtained using the zero temperature analytic form of free electron susceptibility. 

FIG. 7. Half-width at half- maximum of the tunneling conductance peak relative to the spectral 
width at £fc = 0. The inset shows the tunneling conductance at T = 0.05 Tp (dotted line) and 
at T = 0.3 Tp (dashed line) and the tunneling conductance calculated with r e _ e (^,T) replaced 
by r e _ e (0, T). The energy dependence results in a line-shape for the 2D-2D tunneling resonance 
which is not precisely Lorentzian. 
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FIG. 8. Half-width at half-maximum for 2D-2D resonant tunneling peaks: theoretical results 
including local- field corrections to the RPA (full line), theoretical results in the RPA (dashed line), 
Murphy et al. experiment (dotted line). 
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